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ABSTRACT 

The radio galaxy 3C 84 is a representative of 7-ray-bright misaligned active galactic 
nuclei (AGNs) and one of the best laboratories to study the radio properties of the 
sub-pc jet in connection with the 7-ray emission. In order to identify possible radio 
counterparts of the 7-ray emissions in 3C 84, we study the change in structure within 
the central 1 pc and the light curve of sub-pc-size components CI, C2, and C3. We 
search for any correlation between changes in the radio components and the 7-ray 
flares by making use of VLBI and single dish data. Throughout the radio monitor- 
ing spanning over two GeV 7-ray flares detected by the Fermi-LAT and the MAGIC 
Cherenkov Telescope in the periods of 2009 April to May and 2010 June to August, 
total flux density in radio band increases on average. This flux increase mostly origi- 
nates in C3. Although the 7-ray flares span on the timescale of days to weeks, no clear 
correlation with the radio light curve on this timescale is found. Any new prominent 
components and change in morphology associated with the 7-ray flares are not found 
on the VLBI images. 

Key words: galaxies: active- galaxies: jets- galaxies: individual: (NGC 1275, 3C 84). 



1 INTRODUCTION 

With the Large Area Telescope (LAT) on board Fermi 
Gamma-ray Space Telescope, GeV 7-ray emission from 3C 84 
(alias NGC 1275, z=0.0176) was detected for the first time 
(Abdo et al. 2009). It was one of the two misaligned AGNs 
detected by Fermi after four-months operations. More re- 
cently, 7-ray emission is detected in 12 misaligned AGNs 
(Abdo et al. 2010a, Brown & Adams 2012). This gives new 
challenges to the AGN unification scenario. Among these 
misaligned AGNs, 3C 84 exhibited strong variability in its 
7-ray emission. An averaged 7-ray flux during the first four 
months is (2.10±0.23) x 10~ 7 ph cm" 2 s" 1 above 100 MeV. 



* E-mail:hiroshi. nagai@nao.ac.jp 



This 7-ray flux is seven times brighter than the upper limit 
by EGRET. It was claimed that the innermost jet of 3C 84 is 
the most likely source of 7-ray emission because of the time 
variation on the timescales shorter than years to decades 
(Abdo et al. 2009). During the first 2 years, Fermi-LAT de- 
tected 7-ray emission from 3C 84, two episodes of increased 
7-ray activity were observed: one occurred in 2009 Apr-May 
(Kataoka et al. 2010), and the other one occurred in 2010 
June- August (Brown & Adams 2011). In particular, the sec- 
ond flare was detected up to 102.5 GeV and very-high-energy 
7-ray emission with a spectral cut-off around 500 GeV was 
found successively by MAGIC (Aleksic et al. 2011). 

The radio source 3C 84 is known to be a bright ra- 
dio source, and has been studied extensively since the early 
days of radio astronomy. The radio flux has been moni- 
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tored since 1960, and episodes of violent flux increase were 
reported (Kellermann & Pauliny-Toth. 1968; O'Dea et al. 
1984). In mid-1980s, the radio flux became exceptionally 
bright, more than 60 Jy at centimeter wavelengths, and then 
subsequently decreased such that, by the early 2000s, the ra- 
dio flux decreased to ~ 10 Jy (Nagai et al. 2010, hereafter 
Paper I). This flux decrease can be ascribed to the adiabatic 
expansion of radio lobe in the central ~ 5 pc (Asada et al. 
2006). The viewing angle of the jet is estimated to be 11-55° 
to the line of sight (e.g., Asada et al. 2006; Lister et al. 2009). 
Around 2005, it was reported that the radio flux started to 
increase again (Abdo et al. 2009). The VLBI observations 
revealed that this flare originated from within the central 
pc-scale core, accompanying the ejection of a new jet com- 
ponent (Paper I). This new component appeared from the 
south of the core around 2003 (Suzuki et al. 2012, hereafter 
Paper II), and was moving to the position angle ~ 160° 
steadily with slightly changing speed in both parallel and 
perpendicular directions. 

In Paper I, we revealed that the central 1-pc structure 
mainly consists of three bright components CI, C2, and C3. 
In Paper II, we argued that it is difficult to reconcile the ob- 
served GeV 7-ray emission with a one-zoned synchrotron- 
self Compton model in C3 because the observed spectral 
index of C3 between 22 and 43 GHz disagrees with that 
derived from the SED modeling. In order to obtain an addi- 
tional validation to this argument, this paper focuses on the 
comparison between radio light curve and GeV 7-ray flares 
which occurred in the period of 2009 April-May and 2010 
June-August. We argue which component is the most likely 
site of the observed 7-ray emission by comparing radio flux 
variability with the 7-ray flaring events. We are also inter- 
ested in the possible emergence of a new radio component 
associated with the 7-ray flares. 

In this paper, we present the light curve of each compo- 
nent obtained by the Monitoring Of Jets in Active galactic 
nuclei with VLBA Experimental (MOJAVE: Lister et al. 
2009) and detailed structural change obtained by the VLBI 
Exploration Radio Astrometry (VERA). Moreover, the 37- 
GHz data from the Metsahovi have been used to complement 
the MOJAVE light curve. Throughout this paper, we adopt 
H = 70.2 km sec" 1 Mpc" 1 , fl M = 0.27, and Ov = 0.73 (Ko- 
matsu et al. 2011; 1 mas=0.35 pc, 0.1 mas yr" 1 = 0.11c). 



2 OBSERVATION AND DATA REDUCTION 
2.1 The VERA data 

Observations were carried out between 2009 April and 2011 
March for 8 epochs with four VERA stations at 43 GHz. 
Each observation lasted 8 hours, consisting of a number of 
scans in which the length of each scan is 570 seconds, with 
a 30-seconds gap inserted between each scan. Most of the 
scans were assigned for 3C 84 and a few scans were assigned 
for calibration. The total on-source time for 3C 84 is about 7 
hours for each observation. The data were recorded at a rate 
of 128 Mbps for the first 6 epochs and 1024 Mbps for the 



The MOJAVE data archive is maintained at http: / / 
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remaining 2 epochs, with 16-MHz and 128-MHz bandwidth, 
respectively. 

Data reduction was performed using the NRAO Astro- 
nomical Imaging Processing System (AIPS). A standard a 
priori amplitude calibration was performed using the AIPS 
task APCAL based on the measurements of the system tem- 
perature during the observations and the aperture efficiency 
provided by each station for VERA data at 43 GHz. Fringe 
fitting was done using the AIPS task FRING. For the de- 
convolution of synthesized image, we used CLEAN and self- 
calibration technique. Final images were obtained after a 
number of iterations with CLEAN and both phase and am- 
plitude self-calibration using the DIFMAP software package 
(Shepherd 1997). 

After the process described above, we found that the 
VERA 43-GHz was resolved and some portion of flux density 
within the central ~ 1 pc region was missing. Because of this 
missing flux, it was difficult to argue the flux change using 
VERA. We do not use the VERA data for the argument of 
light curve, but they are used to study the structural change 
because the central region is less opaque at 43 GHz than at 
15 GHz, i.e., the frequency of the MOJAVE data. 

2.2 The MOJAVE data 

We imported the calibrated uv-data. of the MOJAVE pro- 
gramme (15 GHz) into the NRAO AIPS package and we 
performed a few phase-only self-calibration iterations and 
flagging some data when the amplitudes were too high or 
too low with respect to the other visibilities. Finally we pro- 
duced high- and low-resolution images in order to study ei- 
ther the core structure or the extended emission. The high- 
resolution image (FWHM=0.55 mas) was produced consid- 
ering only the baselines longer than 50 MA and using pure 
uniform weight, whereas the low-resolution image was ob- 
tained using the baseline shorter than 250 MA and natural 
weight (FWHM=1 mas). 

2.3 The Metsahovi data 

The single dish monitoring data at 37 GHz is adopted from 
the Metsahovi quasar monitoring program. The observations 
were carried out with the Metsahovi 14-m radio telescope. 
A detailed description of the data reduction and analysis is 
given in Terasranta et al. (1998). In order to compare the 
VLBI and single dish light curves, we selected the data from 
2009 January 1 to 2011 November. 



3 MODEL FITS 

In order to investigate if there is any structural change, 
i.e., new component emergence, we deconvolved the VERA 
43-GHz images by using the several components whose 
brightness distribution is elliptical Gaussian. Throughout all 
epochs for VERA data at 43 GHz, overall structures were 
mostly well-represented by three major components CI, C2, 
and C3, which were the same components identified in Pa- 
per I and Paper II. In addition to these three components, 
there was an additional emission bridging between CI and 
C3. We had a difficulty to model this emission: relative flux 
densities and positions of the Gaussian components in this 
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region can be correlated with those of CI and C3, and the 
position and size of CI and C3 were often skewed by these 
components. We operationally fitted this bridging emission 
by using one or two point sources (unresolved components) 
instead of Gaussians (see C4a and C4b in Figure[T]), in order 
to avoid them from such a skewing. The choice of one or two 
point sources was somewhat arbitrary, and it was difficult to 
keep consistency across the epochs. We always tried to main- 
tain consistency of the three main components (CI, C2, and 
C3) even if the position and flux of the point sources C4a 
and C4b did not show smooth change across the epochs. We 
do not use the result from this model fit for the evaluation 
of flux variability because of significant missing flux which is 
mentioned in section [2711 but we use the VERA data for the 
discussion of structural change. Figure [T] shows the CLEAN 
images from the VERA and one example of model- fit image. 
Note that the model-fit image and corresponding CLEAN 
image are essentially similar, demonstrating the validity of 
modeling procedure. 

The analysis of the 15-GHz MOJAVE data was per- 
formed on the image domain instead of on the visibility do- 
main. The flux density, size and position of each component 
was derived by means of the AIPS task JMFIT, which per- 
forms a Gaussian fit to the components on the image plane. 
This approach has been found preferable rather than the 
visibility-based approach due to the presence of extended 
structures. Indeed, the minimum baseline of the VLBA at 15 
GHz (6.5 MA) is sensitive to the extended flux density aris- 
ing from both the northern and southern lobes located about 
15 mas and 10 mas from the core, respectively (e.g., Asada 
et al. 2006). Both components had complex sub-structures 
that make the analysis of the visibility data difficult to carry 
out. On the other hand, the minimum baseline of VERA 
was too long to detect such extended structures, making the 
model fitting of the core components on the visibility data 
more reliable. We use the result from MOJAVE data for the 
discussion of the light curve. 



4 STRUCTURE AND LIGHT CURVE 

Figures [T] and [2] show the total intensity images of MOJAVE 
at 15 GHz and VERA at 43 GHz, respectively. Two bright 
components CI and C3 are separated by about 2 milli-arcsec 
in the north-south direction, while the fainter component C2 
is located on the west of C3 (see Figure[T]). The overall struc- 
ture is very similar to that presented by Paper I and Paper 
II. From the analysis of the data presented in this paper, 
no significant changes in the central structure of 3C 84 has 
been detected. C3 is continuously moving to the south with 
a position angle ~ 170° with respect to CI, as it has already 
been found in previous papers. The positional change of C2 
is somewhat random rather than systematic, but it seems 
that C2 is also moving to the south slightly on average. A 
detailed analysis of the kinematics of the central components 
of 3C 84 will be discussed in a forthcoming paper. 

On average the total flux density increases with time, 
but some rise and fall in flux density on the timescales of a 
month is also seen in the Metsahovi data (Figure [3J. Total 
flux increase shows similar trend with the sum of flux of 
CI, C2, and C3. The most of flux increase originates in C3. 
The increasing level is more than a factor of two throughout 



the monitoring (see Table [TJ . On the other hand, the flux 
variation of C2 and CI is only about 10%. In Paper II it 
was found that C3 showed moderate flux increase until 2008 
and then both CI and C3 showed a rapid flux increase after 
2008. Such a rapid flux increase of CI is not seen anymore. 
Surprisingly, any new component ejection accompanied by 
this flux increase is not detected with the exception of minor 
components C4a and C4b. This is discussed in the section 

El 



5 COMPARISON WITH 7 -RAY EVENTS AND 
RADIO LIGHT CURVE 

In Figure[31 we present the radio data collected by Metsahovi 
and MOJAVE. The pink vertical lines indicate the 7-ray 
flaring periods reported by Kataoka et al. (2010) and Brown 

6 Adams (2011). The first one occurred in the period of 
2009 April-May. The flaring timescale is about 1 month. 
The MOJAVE light curve shows that the flux density of 
CI obtained in 15 days after the peak of this flare (4.98 ± 
0.25 Jy on 2009 May 28) increases slightly as compared to 
that obtained in ~ 4 month before (5.65 ± 0.28 Jy on 2009 
January 30). The flux density of C3 could be also increasing 
between 2009 January 30 (5.56 ± 0.28 Jy) and 2009 May 28 
(5.84 ± 0.29 Jy), but the increasing level is not significant 
as compared to the flux calibration accuracy of VLBA (~ 
5%). On the other hand, no obvious change in radio flux 
density is found from the Metsahovi 37-GHz data on the 
timescale of the 7-ray flaring event. No structural change on 
the VLBI image is also seen before and after the flare despite 
our VERA observations pin down right beginning and end 
of the flare (2009 April 23 and 2009 May 24). 

The second 7-ray flare occurred in the period of 2010 
June- August. Its e- folding rise time is about 1.5 days and a 
subsequent e-folding decay time is about 2.5 days. This flare 
occurred during the increasing activity of total flux density 
at 37 GHz spanning from 2010 to 2011. From Figure [3] it 
is found that this radio flux increase mostly originates in 
C3. However, by focusing on the timescale of days, the light 
curve of total flux does not show any signature of flare. No 
obvious change in VLBI-scale structure before and after the 
flare is seen. 

In the period of both 7-ray flares, no obvious change in 
radio flux density is detected by Metsahovi on the timescale 
of the 7-ray flaring events (days to weeks). If the variability 
in the radio band had been of the same magnitude as in 
the 7-rays, where the flux doubled or tripled its value, the 
single dish observations would have detected such variation. 
Thus, we may conclude that the magnitude of radio variabil- 
ity is smaller than that in 7-ray band or the size of the radio 
counterpart is so small that the radio variability of the flar- 
ing component is mitigated by the constant activity of the 
other components. However, we must note that Metsahovi 
37-GHz light curve shows a small rise and subsequent decay 
in the period from 2009 October to 2009 December (indi- 
cated by the oval of broken line in the small graph in Figure 
3). The timescale of this change is about one month. Right 
after this rise and decay, a small increasing activity which 
occurs on the timescale of one week is detected again (indi- 
cated by the oval of broken dotted line in the small graph in 
Figure 3). Either of these changes in the 37-GHz light curve 
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Table 1. Flux densities of the central components from MOJAVE 
15-GHz data 



Epoch 


CI [Jy] 


C2 [Jy] 


C3 [Jy] 


2009/Jan/30@ 


4.99±0.25 


4.46±0.22 


5.56±0.28 © 


2009/May/28@ 


5.65±0.28 


4.80±0.24 


5.84±0.29 © 


2009/Oct/27@ 


5.78±0.29 


4.74±0.24 


6.30±0.31 @ 


2010/ Aug/6© 


5.80±0.29 


5.30±0.26 


8.61±0.43 @ 


2010/Sep/29@ 


5.51±0.28 


4.24±0.21 


9.35±0.47 © 


2010/Nov/20@ 


5.81±0.29 


4.63±0.23 


10.09±0.50 @ 


2011/Feb/ll@ 


4.84±0.24 


4.22±0.21 


9.85±0.49 © 



Mill ARC SEC 



is possibly related to the 7-ray flare with some time delay, as 
proposed in Leon-Tavares et al. (2011). More recently, an in- 
creasing of radio activity is observed from around 2011 June. 
This probably originates in C3 because the flux density of 
C3 increases between 2011 June 24 and 2011 December 12 
with the significance of 4.5a (Table[T]). This increasing activ- 
ity seems to be difficult to relate to the 7-ray flare because 
the time difference between the two events is very large (~ 1 
year). Further multifrequency analysis to solve this puzzle 
will be presented in a forthcoming paper. 

Where is the radio counterpart of the 7-ray flares? One 
possibility is that the 7-ray flaring region is embedded in the 
self-absorbed core, i.e., CI or further upstream of jet. This 
seems quite natural because the 7-ray time variation is very 
fast, indicating a small emission region size. In particular, 
the resultant emission region size for the case of second flare 
is estimated to be ~ 10~ 3 pc if we assume a Doppler fac- 
tor (S) of 2 (Brown & Adams 2010). The size of ~ 10" 3 pc 
is smaller than the beam size of VLBA at 15 GHz by two 
orders of magnitude. Therefore, it is possible that time vari- 
ation of radio counterpart is mitigated by the contamination 
from other radio emitting region. However, a new jet com- 
ponent associated with the 7-ray flare is likely to be ejected, 
and such a component may be visible on the image some 
time after the 7-ray flare as a result of propagation down to 
optically thin region. If the jet components are ejected from 
the core accompanied by the first and the second 7-ray flares 
with (5 = 2, they should be visible on the image between CI 
and C3 during our monitoring period. As we mentioned in 
section [3] there is a smoothly distributed radio emission be- 
tween CI and C3. This emission can be modeled with two 
point sources. They might be jet components accompanied 
by the 7-ray flares. 

It is worth pointing out for comparison that not all the 
7-ray flare detected so far in blazars have a clear radio coun- 
terpart, as in the case of the GeV-TeV 7-ray flare in 3C 279 
(Abdo et al. 2010b). Another intriguing case in PKS 1510- 
089 is that some 7-ray flares seem to be related to changes 
in radio band (e.g., Marscher et al. 2010; Orienti et al. 2011) 
while others show no relation (e.g., D'Ammando et al. 2009). 
Therefore, the lack of significant changes in radio band for 
3C 84 after the detection of high 7-ray activity leaves the de- 
bate on the region responsible for the high-energy emission 
and its location still open. 



Figure 2. Left: Low-resolution image of MOJAVE 15 GHz on 
2011 December 12. The size of restoring beam is (1 X 1) mas. 
The contours are plotted at the level of (-1, 1, 2, 4, 8, 16, 32, 64, 
128, 256, 512, 1024) X 3.8 mjy/beam, which is the 3-times image 
noise rms. Right: High-resolution image of the central (8 X 8) mas 
region on the same date. The contour levels are same with those 
in low-resolution image. 
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Figure 1. (a)-((g) VERA 43-GHz images and model fit image. The ellipse indicated at the bottom left corner in each image is the 
restoring beam. The size of the restoring beam is (0.63 X 0.37) mas at the position angle of —50°, which is a typical beam size of all 
observing epochs. The contours are plotted at the level of (-1, 1, 2, 4, 8, 16, 32) X 62.6 mjy/beam, which is the 3-times image noise rms 
on 2010 December 27 (f). (h) Model-fit image on 2011 March 10. Symbols indicated by blue lines are the model components obtained 
by the Gaussian model fit. The plus signs are point sources (C4a and C4b) which are placed to model the emission between CI and C3 
(see text in Q. 
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Figure 3. Single dish and VLBI light curves from Metsahovi at 37 GHz and MOJAVE at 15 GHz, respecctively. Red squares, green 
circles, and blue triangles represent the MOJAVE flux densities of CI, C2, and C3, respectively. The sum of flux densities of CI, C2, 
and C3 is indicated by black squares. Single dish light curve by Metsahovi is indicated by cyan diamonds. Two 7-ray flaring events are 
indicated by pink vertical lines. Vertical arrows represent the date of VERA observations. Small graph inserted in overall plot is the 
Metsahovi light curve from 2009.8 to 2010 where the range indicated by the square of thin solid line (see text in 0. Note that Metsahovi 
flux density contains the emission from the entire 3C 84 source. 
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